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ABSTRACT 

The  goal  of  this  thesis  is  to  investigate  the  work  output  of 
a  heat  driven  thermoacoustic  prime  mover  above  onset  of  self- 
oscillation.  The  exponentially  growing  sound  wave,  generated  when 
a  prime  mover  is  initially  "turned  on,"  was  digitally  sampled  for 
a  helium  filled  prime  mover  at  pressures  ranging  from  238  kPa  to 
500  kPa  and  at  temperature  differences  ranging  from  onset  to  400  K. 
This  data  was  then  digitally  filtered  by  a  100  Hz  band  pass  filter 
centered  on  the  prime  mover's  fundamental  frequency.  A  least  mean 
squares  fit  was  applied  to  the  envelope  of  the  filter's  output  in 
order  to  determine  the  temporal  absorption  coefficient  (3.  From  0, 
the  quality  factor  was  computed.  These  quality  factors  were  then 
compared  with  thermoacoustic  theory.  The  agreement  between  the 
theoretical  predictions  and  the  measured  results  is  extremely  good 
at  high  mean  gas  pressures.  As  the  mean  gas  pressure  decreases, 
however,  the  agreement  between  the  theoretical  value  of  onset  and 
the  predicted  slope  of  the  data,  increasingly  deviate. 
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I.   INTRODUCTION 

As  thermoacoustics  is  a  relatively  new  technology,  its 
potential  uses  are  still  being  investigated.  One  such  use, 
which  this  thesis  intends  to  precede,  is  the  use  of  a 
thermoacoustic  prime  mover  as  a  low  freguency  transducer. 
This  possibility  has  been  discussed  by  S.  L.  Garrett  and  T.  B. 
Gabrielson  as  a  viable  alternative  to  the  conventional  sound 
sources  used  today. 

For  conventional  sources,  the  motion  of  a  solid  surface  is 
reguired  to  generate  a  volumetric  velocity.  The  pressure 
amplitude  is  conseguently  limited  by  the  maximum  displacement 
of  the  solid.  The  relationship  between  the  pressure  amplitude 
P  and  the  angular  freguency  of  oscillation  co  for  a  spherical 
source  in  an  infinite  medium  is 

P=Q-£f=Q-pL,  (1) 

2KI  4nr 

where  the  volumetric  velocity  Q  is  given  by 

Q  =   f  vndA.  (2) 

J  surface 

Also,  p  is  the  density  of  the  medium,  a  is  the  speed  of  sound 
in  the  medium,  X  is  the  wavelength  of  the  oscillating  wave,  r 
is  the  range  from  the  source,  and  v,  n,  and  A  are  the  surface 


velocity,  the  unit  vector  normal  to  the  surface,  and  the 
surface  area  of  the  source,  respectively.  Assuming  sinusoidal 
oscillations,  the  surface  velocity  can  be  expressed  as  v  =  wd, 
where  d  is  the  displacement  of  the  surface.  Consequently  Q  is 
proportional  to  a),  the  pressure  amplitude  to  oa2,  and  the 
radiated  acoustic  power  to  gj4.  As  a  result,  as  frequency  is 
decreased,  the  surface  area  must  be  increased  quadratically 
in  order  to  maintain  a  constant  pressure  amplitude  in  the  far 
field  (P  oc  cj2«A)  .  Thermoacoustics  offers  a  solution  to  this 
problem  because  the  driving  force  can  be  applied  directly  to 
the  fluid  instead  of  through  an  oscillating  surface.  With 
only  fluid  motion,  the  limitation  to  the  velocity  no  longer 
comes  from  the  elastic  constraints  of  a  vibrating  solid  but 
from  the  hydrodynamic  (e.g.,  turbulence)  limits  of  the  fluid. 
(Garrett  et  al.,  1990,  pp.  162-163)  The  actual  phenomenon  of 
thermoacoustics  has  been  discussed  at  length  elsewhere.  G.  W. 
Swift's  1988  article  "Thermoacoustic  engines",  for  example, 
gives  a  very  thorough  coverage  of  the  subject.  It  therefore 
will  not  be  explained  in  this  thesis. 

Thermoacoustics  research  at  the  Naval  Postgraduate  School 
has  been  directed  at  thoroughly  verifying  thermoacoustic 
theory.  In  1989,  Hsiao-Tseng  Lin  studied  a  thermoacoustic 
prime  mover  below  onset  (Lin,  1989) .  The  conclusions  of  his 
thesis  indicated  that  theory  based  on  short-stack  and  boundary 
layer  approximations  did  not  hold  valid  in  all  cases.  He  also 


found  that  above  onset,  the  acoustic  amplitude  quickly  became 
quite  large,  and  the  waveforms  highly  nonlinear. 

The  next  logical  progression  in  this  research  is  to  study 
a  prime  mover  above  onset  while  the  acoustic  amplitude  is 
small  and  to  apply  a  non  short  stack,  boundary  layer  theory. 
It  is  intended  that  the  information  gained  will  help  lead  to 
an  understanding  of  prime  movers  at  large  amplitudes.  The 
purpose  of  this  thesis  is  to  take  the  first  step  by  looking  at 
the  initial  growth  of  the  oscillations  as  they  build  up  to 
steady  state. 


II.   THEORY 

A.   INTRODUCTION 

A  prime  mover  can  be  divided  into  five  major  regions. 
These  five  regions,  as  shown  in  Figure  1,  are  the  ambient  end, 
the  ambient  heat  exchanger,  the  stack,  the  hot  heat  exchanger, 
and  the  hot  end.  In  typical  configurations,  the  ambient  and 
hot  ends  are  empty  capped  tubes  maintained  at  constant 
temperature.  Since  we  chose  to  apply  a  differential 
temperature  to  the  prime  mover  by  heating  one  end  while 
keeping  the  other  end  at  ambient  temperature,  hot  and  ambient 
are  used  to  describe  the  different  regions.  The  ambient  and 
hot  heat  exchangers  both  consist  of  a  series  of  plates  used  to 
maintain  the  opposing  ends  of  the  stack  region  at  their 
respective  temperatures.  Finally,  the  stack  region  is  where 
the  conversion  of  thermal  to  acoustic  energy  takes  place. 

The  purpose  of  this  thesis  is  to  compare  the  measured  work 
output  of  a  prime  mover  above  onset  of  self-oscillation  to  the 
theoretical  value.  Unfortunately,  it  is  not  a  simple  matter 
to  directly  measure  the  work  output.  What  is  easily  measured 
though,  is  the  acoustic  signal  generated  by  the  prime  mover. 
The  evolution  of  this  signal  can  be  related  to  the  work  output 
as  discussed  below. 
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Figure  1  Major  regions  and  components  of  a  prime  mover 


When  a  resonator  is  excited  and  then  allowed  to  decay,  the 
amplitude  of  the  oscillations  decay  according  to 

(3) 

g-ut/2Pgjot  *      ' 


where  t  is  time,  Q  is  the  quality  factor  of  the  resonator  and 
j2  =  -1.  (Swift,  1988,  p.  1162)  In  terms  of  the  temporal 
absorption  coefficient  0,    this  exponential  behavior  is 

e-9tej»ti  (4) 

where  (3  =  co/2Q.  The  initial  buildup  of  oscillations  in  a 
prime  mover  can  be  characterized  by  a  negative  value  of  /3 . 
Thus,  the  amplitude  initially  grows  exponentially.  As  the 
oscillations  become  larger  in  amplitude,  however,  the  value  of 
the  net  absorption  coefficient  changes  and  eventually  equals 
zero.  This  results  in  steady  state  oscillations  of  finite 
amplitude  as  can  be  seen  in  Figure  2.  It  is  the  initial 
growth  of  the  signal  that  is  of  interest  to  us. 

|3,  which  is  equivalent  to  the  inverse  of  the  rise  time, 
can  be  directly  obtained  from  analyzing  the  envelope  of  the 
acoustic  signal  generated  by  the  prime  mover.  Along  with  /3, 
the  frequency  of  oscillation  is  also  measured  so  that  1/Q  may 
be  computed.  Recall  that  1/Q  =  2/3/w  =  &/nf,  where  f  is  the 
oscillation  frequency  measured  in  hertz.  The  choice  of  using 
1/Q  as  the  parameter  of  interest  comes  from  the  fact  that  it 
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Figure  2  Signal  which  shows  the  changing  of  0  from  an 
initial  negative  value  (evidenced  by  the  exponential 
buildup)  ,  to  zero  (as  indicated  by  the  steady  state 
result) . 


is  directly  proportional  to  the  work  output  of  the  prime 
mover.  As  discussed  by  Swift,  1/Q  relates  to  the  average 
dissipated  acoustic  power  E  of  a  prime  mover  by 

1  -  *  (S, 


Q         uE 


ST 


where  EST  is  the  energy  stored  in  the  resonator  (Swift,  1988, 
p.  1162) .  Rather  than  using  the  dissipated  acoustic  power  E, 
we  prefer  to  work  with  the  generated  acoustic  power  W.  W  is 
commonly  referred  to  as  the  work  output.  By  convention,  W  is 
positive  when  power  is  generated  while  E  is  positive  when 
power  is  dissipated.   Therefore,  E  =  -W,  and 

1  -   "*  (6) 


Q         (j)E 


ST 


B.   AVERAGE  GENERATED  ACOUSTIC  POWER,  WsUck 

The  acoustic  power  generated  in  the  stack  section  of  a 
prime  mover  cannot  flow  or  escape  in  a  direction  perpendicular 
to  the  stack  plates,  the  y-direction.  Energy  transport  is 
allowed  however,  in  the  x-direction  which  is  parallel  to  both 
the  plates  and  the  temperature  gradient.  A  measurement  of  the 
difference  in  average  acoustic  intensity  between  the  two  ends 
of  the  stack  can  therefore  be  used  to  determine  the  acoustic 
power  generated  by  the  stack. 


The  acoustic  power  generated  by  a  short  segment  of  stack 
is  given  by 

Ktack    =    S[(pu)     right      ~     {pU)      left     ] 
end  of  end  of 

segment  segment  .  _ . 

=    SAX-p-  (PU) 

dx 


where  S  is  the  cross  sectional  area  of  the  fluid  in  the  stack, 
pu  is  the  time  average  acoustic  intensity,  and  Ax  is  the 
length  of  the  short  segment.  Note  that  SAx  is  the  total 
volume  of  the  fluid  in  the  stack  region.  By  dividing  through 
by  this  volume,  the  power  per  unit  volume  is  obtained; 


W. 

w  = 


"?  ■  4.  <P5>  -  ^elP^  ♦  ^  M  ■     <8> 


5ax    dx  2  1     dx  dx 

where  a  tilde  denotes  complex  conjugate,  and  p,  and  u, 
represent  the  spatial  dependence  of  the  pressure  and  velocity, 
respectively,  (i.e.  p(x,t)  =  p^"*  and  u(x,t)  =  u^"')  .  Also, 
the  angled  brackets  denote  the  spatial  average  of  velocity 
across  the  fluid  gap  between  the  plates.  It  is  necessary  to 
spatially  average  the  acoustic  velocity  across  the  fluid  gap 
because,  with  the  presence  of  viscosity,  u  is  also  y- 
dependent.  Expressions  for  the  four  terms  on  the  right  hand 
side  of  this  equation  will  be  derived  next  in  order  to  arrive 
at  a  complete  equation  for  work  output.  (Swift,  1988,  p.  1159) 


Since  the  prime  mover  is  essentially  a  pipe  with  rigid 
ends,  the  pressure  behaves  sinusoidally  like  a  cosine  wave. 

Pj_  =  PAcos  (kx)  ,  (9) 

where  PA  is  the  peak  pressure  amplitude,  k  is  the  wave  number, 
and  x  =  0  at  the  capped  end  of  the  ambient  end. 

Similarly,  it  can  be  shown  that  the  velocity  in  the  stack 
is  represented  by 

Uj.  =  -j'-^-(l  +  —  )sin(Jcx)  s  jUla,  (10) 

Pm«3      Vo 


where  pm  is  the  mean  density  of  the  fluid,  1  is  one-half  the 
plate  thickness,  and  y0  is  one-half  the  plate  spacing.  (Swift, 
1988,  p.  1157) 

From  the  equation  of  motion  for  a  compressible,  viscous 
fluid,  Swift  has  shown  that 


dPi   =  -J(Jpm(  Ui  )  =  upm(  u±_  )  (11) 


dx  1  •■  fv  1  -  fv    ' 


where 


tanhld  *  j)yjij 

(1  *  j)yJK  (12' 


■\> 


10 


and  6„  =  v  2/i/poj  is  the  fluid's  viscous  penetration  depth  and 
ju  is  the  dynamic  viscosity  of  the  fluid  (Swift,  1988,  p. 
1159)  . 

Finally,  the  wave  equation  can  now  be  used  to  solve  for 
d<u,>/dx.  (Swift,  1988,  p.  1157) 


1  +  €.    x    a)2  dx         pm       dx 


'  s 


2 


f,  -  f*         vr,  dPx  _ 


(13) 


a 

a)2  (1  ■■  o)  (1  +  €s)  Tm     dx 


=   0 


where 


.  =  P^p5Ktanh[(l  *  j)y0/6K] 
pscs6stanh[(l  +  j)l/b3]   ' 


cp  and  cs  are  the  isobaric  specific  heats  per  unit  mass  for  the 
fluid  and  the  solid  (i.e.  the  plate),  8K  =  V  2K/pmcpw  and  <5S  = 
vaKj/PjCjCj  are  the  thermal  penetration  depths  in  the  fluid  and 
the  solid,  ps  is  the  density  of  the  solid,  and  K  and  K,  are  the 
thermal  conductivity  of  the  fluid  and  the  solid.  Also,  7  is 
ratio  of  the  fluid's  isobaric  to  isochoric  specific  heats,  a 
is  the  fluid's  Prandtl  number,  VTm  is  the  mean-temperature 
gradient  in  the  x-direction,  Tm  is  the  mean  temperature,  and 
fK  is  a  function  similar  to  f„.  Substituting  Eq.  (11)  into  the 
second  term  of  the  wave  equation  yields  the  following; 
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cK  ut  ) 


dx 


=  J 


fK    -    ^v 


VTm  dPl 


"Pm 


(l  -  a)  (1  +  €.) 


dx 


-  J 


0) 


Pm<3 


<-4^>p, 


(15) 


Applying  these  four  expressions,  Eqs.  (9) ,  (10)  ,  (11) ,  and 
(15),  to  Eq.  (8)  gives  the  equation  we  have  been  working  for; 


w  = 


^4 
2 


Pfo 


P»a 


[J( 


f.-£. 


l  + 


yc 


Vr 


5--/ r-7 __2?tan(icx)  +  1 

1  -  £„  (1  -  0)  (1  +  €s)  kTm 


(16) 


1    V2 


(1  +  — ) 

+  -1 -2L)  cos2  (kx)    -  j 


1  +  €. 


1  -  fL 


sin2  (icxr) 


Now  with  w  known,  computing  the  total  acoustic  power 
simply  requires  integrating  over  the  volume  of  the  fluid 
within  the  stack; 


"scac*   =  f*dV   "  f*SdX. 


(17) 


C.   TOTAL  ACOUSTIC  POWER,  W 

To  compute  1/Q,  the  total  acoustic  power  W  is  needed.  In 
the  previous  section,  only  an  expression  for  the  work  done  by 
the  prime  mover's  stack  was  derived.    The  four  remaining 
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sections  of  the  prime  mover  also  need  to  enter  into  the 
equation.  The  total  work  done  by  the  prime  mover  can  be 
divided  into  the  sum  of  the  work  done  by  or  on  each  of  these 
five  individual  sections, 

W  ~    ^amb  end    +    ^ambhx   +    ^ stack    +    ^hot  hx    +    ^ hot  end'  (18) 

It  should  be  noted  that,  with  the  exception  of  W,uck,  these 
terms  represent  losses.  WsUck  can  represent  either  loss  or  gain 
depending  on  the  value  of  the  temperature  gradient. 

The  equation  for  w  was  derived  for  the  stack  region  but  it 
actually  applies  to  each  of  the  remaining  four  regions.  The 
parameters  in  w  simply  vary  depending  on  what  section  of  the 
prime  mover  x  happens  to  be  in.  In  all  four  remaining  regions 
for  example,  VTm  =  0  since  there  is  no  temperature  gradient 
anywhere  except  in  the  stack  region.  Likewise,  1  =  0  in  the 
ambient  and  hot  ends  since  these  regions  do  not  contain  any 
plates  and  the  functions  f,  and  f,  take  on  slightly  different 
forms  as  well. 

The  total  acoustic  power  is  consequently  one  integral  from 
the  beginning  of  the  ambient  end  to  the  end  of  the  hot  end, 


W  =   f L  w(x)S(x)dx.  (19) 

Jjf-0 


Within  this  integral  numerous  parameters  vary  with  respect  to 
x  depending  upon  which  region  of  the  prime  mover  x  happens  to 
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be  in  (e.g.,  perimeter  size,  plate  thickness,  plate  spacing). 
Also,  within  the  given  region  of  the  stack  there  are  several 
parameters  which  vary  as  a  function  of  x  because  of  the 
temperature  change  along  the  stack  (e.g.,  sound  speed,  thermal 
conductivity  of  the  fluid,  viscosity,  etc.).  All  such 
parameters  are  treated  as  functions  of  x  so  that  a  single 
integration  yields  the  total  work  done  by  the  prime  mover. 

D.   STORED  ENERGY,  EST 

The  final  quantity  needed  to  calculate  1/Q  is  the  energy 
stored  in  the  resonator.  This  is  found  by  integrating  the 
time  averaged  energy  density  over  the  volume  of  the  resonator. 


EST   =  fedV.  (20) 


where  the  instantaneous  energy  density  is  given  by 


5  =  h.i-£z +  u2>< 


2  pW 


(21) 


and  its  time  average  by 


♦-ip.<UiV.  (22) 


Pm<3 


4.'« 
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Again,  as  in  the  case  of  acoustic  power,  the  total  stored 

energy  can  be  represented  as  the  sum  of  the  stored  energy  in 
each  of  the  prime  mover's  five  regions, 

•3J                  i-'a«rf>Md                  i>Iambhx                 * 1  stack                 iJ4otJw  ^ 1  hot  and 


Choosing  appropriate  expressions  for  p,  and  u,'  and  integrating 
the  time  averaged  energy  density  over  each  of  the  five  regions 
produces: 

fioT  ~  7 1     X 


'ST 


A  p  I       amb  end 

4   n    a21 


Pma 


+   r         [Co3>(kx  +  V  +   (1  +  _l)sin2(^x  +  4))]dx 

Ja^bhx          1  +    1  yQ 

yQ 

+   r         [cos^jcx^)  +   (1  +  ±)sin2{kx  +  4>)]dx 

J  stack              1    +    A  -Vo 

+   r         [co3Mjcx_L»)  +   (1  +  A)sin2(icx.4))]dx 

J  hot  hx            -           1  y„ 


xAot  endj 


where  <p  is  the  phase  of  the  standing  wave  at  the  entrance  of 
the  region.  Finally,  combining  Eqs.  (16),  (19)  and  (24)  with 
Eq.  (6)  yields  the  desired  parameter,  1/Q. 
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III.   EXPERIMENTAL  SETUP 

A.   THE  PRIME  MOVER  TUBE 

The  physical  construction  of  the  prime  mover  assembly  was 
previously  shown  in  Figure  1.  Its  internal  diameter  is  3.82 
cm  and  its  overall  internal  length  is  1.0005  m.  A  complete 
listing  of  all  important  prime  mover  measurements  is  found  in 
Appendix  A. 

1.  The  Hot  Heat  Exchanger  and  the  Hot  End 

The  hot  heat  exchanger  and  the  hot  end  tube  are  both 
constructed  of  nickel  as  is  the  end-cap  of  the  hot  end  tube. 
This  end-cap  is  welded  to  the  hot  end  tube  and  presents  a 
rigid  boundary  for  sound  waves  generated  by  the  prime  mover. 
Through  the  center  of  this  end-cap  a  small  hole  has  been 
drilled  to  accommodate  a  thermocouple  probe  for  sensing  the 
hot  heat  exchanger  temperature.  To  prevent  leakage  of  gas 
through  this  hole,  the  thermocouple  probe  is  epoxied  in  place. 

The  hot  end  tube  is  0.63  cm  thick  and  has  an  internal 
length  of  5.50  cm.  Its  internal  diameter,  as  are  all  of  the 
individual  components  which  make  up  the  prime  mover,  is  3.8  2 
cm.  The  hot  heat  exchanger  assembly  consists  of  25,  0.04  5  cm 
thick  and  0.762  cm  long  nickel  plates.  The  spacing  between 
each  plate  is  0.104  cm.  To  minimize  deformation  from  high 
temperatures,  spacers  made  from  3  04  stainless  steel  with  a 
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diameter  of  0.031  cm  are  welded  to  one  side  of  each  plate. 
The  hot  heat  exchanger  assembly  is  mounted  to  the  open  end  of 
the  hot  end  tube. 

2 .  The  Stack 

The  prime  mover  stack  is  housed  in  a  304  stainless 
steel  tube  of  0.05  cm  thickness  with  a  length  of  3.50  cm.  It 
is  comprised  of  35,  0.028  cm  thick  304  stainless  steel  plates 
and  is  the  same  length  as  the  tube.  The  plate  spacing  in  the 
stack  region  is  0.077  cm  and,  as  in  the  hot  heat  exchanger, 
spacers  are  attached  to  each  plate  to  minimize  deformation. 
The  stack  and  the  hot  heat  exchanger  are  assembled  such  that 
their  plates  are  at  near  right  angles  to  one  another  so  that 
the  plates  of  the  heat  exchanger  do  not  coincide  with  the  gaps 
in  the  stack.  Such  a  misalignment  would  result  in  the 
blockage  of  fluid  flow. 

3.  The  Ambient  Heat  Exchanger  and  Ambient  End 

The  ambient  heat  exchanger  consists  of  two  identical 
copper  stacks  each  containing  25,  1.02  cm  long,  0.045  cm  thick 
plates.  The  gap  between  each  plate  is  0.104  cm  and  the 
spacing  between  the  two  individual  stacks  is  0.15  cm.  These 
two  stacks  are  housed  within  a  brass  water  jacket  and  the  hot 
end/stack  assembly  is  soldered  to  one  end  of  the  jacket.  The 
other  end  of  the  jacket  is  soldered  to  the  ambient  end  of  the 
prime  mover  which  is  constructed  of  an  88.11  cm  long  copper 
tube.  A  butterfly  valve  is  positioned  within  this  tube,  38  cm 
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from  the  ambient  heat  exchanger.  The  purpose  of  this  valve  is 
to  shut  off  the  sound  wave  generated  by  the  prime  mover.  By 
closing  the  valve  the  internal  length  of  the  prime  mover  is 
effectively  shortened  reguiring  a  much  higher  differential 
temperature  to  drive  the  prime  mover  above  onset.  In  this 
way,  a  differential  temperature  which  places  the  prime  mover 
above  onset  can  be  established  without  sound  being  generated. 
By  opening  the  valve  the  prime  mover  is  instantaneously  placed 
above  onset  (turned  on)  and  the  initial  pressure  wave  build  up 
may  be  observed  and  measured.  Lastly,  the  ambient  end  is 
flanged  and  a  brass  end-cap  which  houses  a  pressure  sensing 
transducer  is  fastened  to  it. 

B.   CONTROL  AND  MEASUREMENT  EQUIPMENT 
1.  Temperature 

The  temperature  of  the  hot  heat  exchanger  is 
controlled  by  an  Omega  Engineering  HBA  Model  202040  heater 
powered  by  a  GR  Type  100-Q  variable  transformer.  The  heater 
coils  are  tightly  wrapped  and  cover  the  entire  hot  end  and  hot 
heat  exchanger  assembly  to  provide  a  constant  "hot" 
temperature  to  this  portion  of  the  prime  mover.  As  previously 
mentioned,  a  type  K  thermocouple  entering  through  the  hot  end 
end-cap  is  positioned  in  contact  with  the  hot  heat  exchanger 
to  measure  the  hot  end  temperature.  During  operation  the  hot 
end  assembly  is  surrounded  with  insulation  to  minimize  heat 
loss. 
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Ambient  end  temperature  is  controlled  through  the  use 
of  a  Neslab  RTE-110  constant  temperature  bath.  The  bath  is 
used  to  cycle  constant  temperature  water  through  the  water 
jacket  surrounding  the  ambient  heat  exchanger  and  then  through 
plastic  tubing  which  is  wrapped  around  the  entire  length  of 
the  ambient  end.  Temperature  is  monitored  via  three  type  E 
thermocouples  which  are  glued  to  the  top,  middle,  and  bottom 
sections  of  the  ambient  end.  This  method  of  ambient 
temperature  control  maintains  the  ambient  end  temperatures 
within  a  1°C  spread.  The  hot  end  thermocouple  along  with  the 
three  ambient  end  thermocouples  are  monitored  using  a  Keithley 
74  0  System  scanning  thermometer. 
2.  Pressure 

Mean  pressure  inside  the  prime  mover  is  sensed  via  a 
connection  which  penetrates  the  prime  mover  between  the  two 
ambient  heat  exchanger  stacks.  This  connection  is  valved  so 
that  the  pressure  sensing  line  may  be  isolated  from  the  prime 
mover  when  oscillations  are  being  measured.  Leaving  the  valve 
open  increases  the  losses  of  the  prime  mover  and  actually  puts 
it  below  onset  for  the  differential  temperatures  examined. 
Pressure  is  measured  from  this  sensing  line  by  an  Omega  PX3  04- 
3  00AV  pressure  transducer  and  an  HP-3478A  multimeter  is  used 
to  display  the  voltage  output  of  the  transducer.  Pressurizing 
the  prime  mover  is  done  by  first  evacuating  the  prime  mover 
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chamber  with  a  vacuum  pump  and  then  filling  it  with  helium 
gas.   All  this  is  done  via  the  pressure  sensing  line. 
3.  Signal 

The  acoustic  signal  generated  by  a  differential 
temperature  applied  to  the  prime  mover  is  sensed  by  an  ENDEVCO 
Model  8510B-5  piezoresistive  pressure  transducer  housed  within 
the  ambient  end  end-cap.  The  linear  response  range  of  the 
transducer  is  ±5  psi.  A  high  impedance  leak  is  provided 
between  the  resonator  and  the  back  volume  of  the  transducer  to 
eliminate  the  dc  pressure  difference  on  the  transducer  while 
having  negligible  effect  on  the  acoustic  pressure  differences. 
The  leak  is  a  0.004"  internal  diameter,  approximately  10  cm 
long  piece  of  copper-nickel  capillary  tubing.  The  signal  from 
the  transducer  is  routed  through  a  locally  constructed 
differential  amplifier  which  boosts  the  signal  by  a  factor  of 
ten.  This  amplified  signal  is  then  applied  as  input  to  a 
Nicolet  310  digital  oscilloscope  for  data  collection  and 
storage. 
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IV.   EXPERIMENTAL  PROCEDURE  AND  SIGNAL  PROCESSING 

A.   DATA  COLLECTION 

The  first  step  in  obtaining  prime  mover  response  data  was 
to  evacuate  the  resonator.  This  was  accomplished  by  pumping 
down  the  prime  mover  via  its  pressure  sensing  line  for  several 
hours.  Once  fully  evacuated,  the  vacuum  pump  was  turned  off 
and  isolated  from  the  prime  mover.  The  tube's  internal 
pressure  was  then  recorded  and  this  value  was  used  as  the 
calibration  value  for  the  pressure  sensor.  All  future 
pressure  readings  were  corrected  by  this  amount.  The  prime 
mover  was  then  filled  with  the  desired  pressure  of  helium  to 
within  a  tolerance  of  ±  1  kPa. 

For  temperature  control,  the  constant  temperature  bath  was 
started  and  the  variable  transformer  was  turned  up  to  power 
the  heater  element.  Once  a  steady  state  differential 
temperature  of  interest  was  reached,  the  prime  mover  was  ready 
for  data  collection. 

As  mentioned  previously,  a  Nicolet  310  digital 
oscilloscope  (NIC-310)  was  used  to  sample  the  waveform 
generated  by  the  prime  mover.  Since  the  data  was  being 
digitally  sampled  it  was  important  to  at  least  satisfy  the 
Nyquist  sampling  rate  of  2fmix  in  order  to  prevent  aliasing, 
where  f^  is  the  highest  frequency  of  interest  (Strumm,  1988, 
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p.  54) .  For  the  prime  mover,  we  are  only  concerned  with  the 
rise  time  of  the  fundamental  component  of  the  acoustic  signal 
which  has  a  frequency  of  approximately  500  Hz.  Consequently 
we  must  sample  at  a  minimum  of  1000  Hz.  To  be  conservative  a 
sampling  frequency  of  four  times  the  fundamental  was  chosen  as 
the  minimum  acceptable  sampling  frequency.  As  it  turns  out, 
the  closest  sampling  frequency  offered  by  the  NIC-310  without 
violating  this  threshold  is  5000  Hz  or  200  /xs  per  data  point. 
With  the  4000  point  data  capacity  of  the  NIC-310,  this  process 
enabled  the  recording  of  sample  records  of  up  to  0.8  seconds 
in  duration. 

Data  collection  was  done  by  placing  the  onset  control 
valve  in  the  closed  position  to  turn  off  the  sound  generated 
by  the  tube.  After  thermal  equilibrium  was  achieved,  the 
valve  was  opened.  When  a  signal  could  be  seen  over  the 
background  noise  present,  the  oscilloscope  record  was  frozen 
and  saved  on  a  3.5"  floppy  disk  for  later  analysis.  Several 
such  records  were  taken  for  each  pressure  and  temperature 
combination  of  interest. 

B.   DIGITAL  SIGNAL  PROCESSING 

Analysis  of  the  accumulated  data  was  done  with  MathWorks' 
program  MATLAB.  The  actual  data  records  stored  by  the  NIC- 
310,  however,  were  in  binary  code  so,  to  make  the  data  records 
ready  for  processing,  they  were  simply  converted  into  ASCII 
files. 
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Initially,  for  each  data  record,  a  Discrete  Fourier 
Transform  (DFT)  was  computed  and  the  fundamental  frequency  of 
the  record  was  identified.  This  frequency  was  then  used  to 
translate  the  frequency  spectrum  of  the  signal  so  that  the 
fundamental  frequency  was  now  located  in  the  first  bin  of  the 
DFT.  This  was  accomplished  by  multiplying  the  signal  by  a 
complex  exponential  in  the  following  fashion; 

xnew(n)    =  e3{2n/mMn  xold{n)  ,  (25) 

where  N  is  the  number  of  points  contained  in  the  data  record, 
n  is  an  integer  denoting  the  position  within  a  record,  and  M 
is  the  number  of  bins  which  every  frequency  bin  in  the  DFT  of 
xoW(n)  is  to  be  circularly  shifted.  For  example,  if  the 
fundamental  frequency  of  the  DFT  of  xold(n)  resides  in  bin  M  + 
1,  then  a  DFT  of  xncw(n)  will  have  the  fundamental  frequency 
located  in  its  first  bin.  (Strumm,  1988,  p.  429)  The  DFT  of 
a  typical  data  set  is  displayed  in  Figure  3.  Note  that  there 
are  two  peaks  in  this  DFT.  This  is  because  it  is  the  DFT  of 
a  real  (non  complex)  signal.  Since  real  signals  are 
sinusoidal,  (i.e.  the  sum  of  two  exponential  functions  of 
equal  magnitude,  but  with  frequencies  of  opposite  sign) ,  two 
fundamental  frequencies  result.  One  represents  the  positive 
frequency  of  the  sine  wave,  the  first  half  of  the  DFT,  and  the 
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Figure  3   Discrete  Fourier  Transform  (DFT)  of  a  typical 
data  record  showing  the  presence  of  a  500  Hz  signal. 
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other  represents  the  negative  frequency  of  the  sine  wave,  the 
second  half  of  the  DFT. 

Now,  with  the  fundamental  frequency  effectively  translated 
to  dc  level  (first  bin),  a  low  pass  filter  can  be  applied  to 
isolate  the  desired  signal.  The  low  pass  digital  filter 
chosen  for  this  task  was  a  fifth  order  Butterworth  filter  with 
a  bandwidth  of  50  Hz.  This  procedure  is  essentially  the  same 
as  filtering  the  original  data  with  a  100  Hz  band  pass  filter 
centered  at  the  fundamental  frequency  of  the  data  record.  By 
translating  the  data  first,  however,  we  immediately  arrive  at 
the  envelope  of  the  waveform  since  its  frequency  has  been 
"converted"  to  dc.  Figures  4  and  5  show  a  typical  data  record 
before  and  after  processing,  respectively.  Note  that 
processing  the  signal  in  this  manner  has  the  effect  of  halving 
the  magnitude  of  the  original  signal. 

In  most  cases  the  4000  point  DFT  revealed  that  only  the 
fundamental  frequency  was  present.  As  the  prime  mover  was 
driven  higher  (-50K)  above  onset,  however,  the  second  and  even 
the  third  harmonics  were  observed.  To  ensure  that  only  the 
affects  caused  by  the  fundamental  were  measured,  DFT's  were 
taken  of  the  earliest  parts  of  these  data  records  to  eliminate 
the  second  harmonic.  This  is  important  because,  as  the  second 
and  higher  harmonics  begin  to  grow,  they  take  energy  away  from 
the  first  harmonic  and  cause  /3  to  change.  Such  losses  are  not 
accounted  for  by  our  theory.   This  effect  can  be  seen  in 
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Figure  4    Typical  data  record  before  processing 
(Sampled  at  5000  samples  per  second.) 
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Figures  6  and  7  where  the  first  2000  data  points  contain  only 
the  fundamental  frequency  while  the  entire  4000  point  sample 
shows  evidence  of  the  first  three  harmonics.  Note  that  for 
the  first  0.4  seconds  (2000  samples),  the  slope  of  the 
logarithm  of  the  envelope  is  essentially  constant  (Figure  7). 
This  slope  is  equivalent  to  (3  since  d/3/dt  is  approximately 
zero.  As  the  harmonics  begin  to  appear,  however,  the 
magnitude  of  0  noticeably  decreases  as  indicated  by  the 
lessening  of  the  slope  on  the  logarithmic  plot  of  the  data. 
The  2000  and  4000  point  DFT's  of  this  data  are  shown  in 
Figures  8  and  9. 

As  with  all  digital  filtering,  the  end  conditions  are  the 
most  difficult  to  match.  Because  of  this  the  initial  and 
final  points  of  a  filtered  signal  need  to  be  discarded.  To 
determine  the  actual  range  of  the  filtered  data  points  which 
are  "influenced"  by  the  filtering  process,  test  signals  were 
generated  of  the  following  form; 

TestSignal  =   0.001  ept  sin(27i  500t)  .        <26) 

The  use  of  a  scaling  factor  of  0.001  and  a  frequency  of  500  Hz 
were  made  to  generate  signals  which  were  characteristic  of  the 
data  being  analyzed.  Likewise,  varying  values  of  (3  were  used 
to  cover  the  range  of  0's  that  were  being  measured.  The 
results  show  that  smaller  values  of  (3  produce  larger  relative 
fluctuations  in  the  filtered  waveform.  Figures  10  and  11  show 
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Figure  6  Filtered  signal  (solid  line)  which 
experiences  the  growth  of  higher  harmonics  after 
the  first  0.4  seconds.  The  dashed  line  is  a  least 
mean  squares  fit  to  the  first  0.4  seconds  of  data. 
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Figure  7  Logarithmic  plot  of  signal  shown  in 
Figure  6  (solid  line) .  The  dashed  line  is  a  least 
mean  squares  fit  to  the  first  0.4  seconds  of  data. 
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Figure  8  A  2000  point  DFT  of  the  signal  shown  in 
Figures  6  and  7,  showing  that  only  the  fundamental 
frequency  is  present  in  the  first  0.4  seconds  of 
the  signal. 
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Figure  9  A  4000  point  DFT  of  the  signal  shown  in 
Figures  6  and  7 ,  showing  that  three  harmonics  are 
present  in  the  entire  0.8  second  sample. 
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Figure  10  A  filtered  test  signal  with  /3  =  0.5 
showing  the  endpoint  influence  caused  by  the 
filter. 
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Figure  11  A  logarithmic  plot  of  Figure  10  showing 
the  endpoint  influence  of  the  filter  on  a  smooth 
test  signal. 
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both  the  linear  and  logarithmic  presentations  for  a  test 
signal  ((3  =  0.5)  that  was  processed  in  accordance  with  the 
above  procedure.  As  these  plots  indicate,  the  filter's 
endpoint  effect  extends  to  the  first  and  last  several  hundred 
data  points.  To  verify  this,  a  least  mean  sguares  fit  was 
applied  to  the  natural  logarithm  of  the  test  signal  shown.  In 
an  attempt  to  eliminate  the  filter's  endpoint  effect  the  first 
and  last  300  points  were  discarded.  The  result  was  that  (3 
could  be  determined  to  within  one  millionth  of  one  percent,  a 
negligible  error.  For  data  records  which  contained  only  the 
fundamental  harmonic,  this  criteria  of  discarding  the  initial 
and  final  300  data  points,  after  filtering,  was  used  to 
prevent  the  filter  from  adversely  affecting  the  processed 
results. 

In  the  case  where  a  signal  was  observed  to  contain 
multiple  harmonics,  the  filtered  data  was  truncated  at  some 
value  less  than  3700.  The  actual  value  was  determined  by  two 
criteria.  First,  as  previously  mentioned,  the  DFT  was  used  to 
determine  at  what  point  in  a  data  record,  the  buildup  of  the 
second  harmonic  began.  This  value  marks  the  maximum 
acceptable  truncation  level.  A  least  mean  squares  fit  was 
then  applied  to  the  logarithm  of  the  filtered  signal  truncated 
at  this  value.  To  evaluate  the  quality  of  this  truncation, 
the  filtered  signal  was  plotted  along  with  the  least  mean 
squares  fit  conducted  over  this  truncated  range.  If  the  two 
curves  did  not  visibly  agree  over  the  entire  truncated  range 
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then  it  was  obvious  that  the  fundamental  was  not  completely 
isolated.  This  can  occur  because  the  initial  buildup  of  a 
harmonic  is  difficult  to  detect.  Because  we  are  dealing  with 
very  narrow  band  frequencies,  they  will  initially  go 
undetected  since  the  frequency  bins  are  1.25  hertz  wide  for  a 
0.8  second  sample.  Also,  as  less  data  is  used  in  performing 
a  DFT  in  an  attempt  to  isolate  the  fundamental  frequency,  the 
frequency  bins  become  wider.  In  cases  where  the  results  of 
the  initial  truncation  were  unsatisfactory,  a  smaller  portion 
of  the  sample  was  chosen  and  the  data  was  reprocessed  under 
this  same  criteria  until  an  agreement  was  evident.  Appendix 
B  gives  a  tabulation  of  each  data  record  analyzed  and  the 
range  of  the  points  that  were  used  for  determining  /?. 

This  chosen  segment  of  the  filtered  data  record  now 
represents  the  exponential  behavior  of  the  fundamental 
frequency  alone  and,  assuming  that  d/3/dt  is  small,  the  natural 
log  of  this  data  is  a  line  whose  slope  is  equivalent  to  0.  By 
performing  a  least  mean  squares  fit  to  this  data,  a  value  for 
0  can  be  obtained.  Figures  12  through  15  show  four  typical 
data  records.  They  reveal  that  jS  does  not  noticeably  change 
over  the  time  that  the  data  is  observed  when  only  the  first 
harmonic  is  present.  The  uncertainty  in  the  slope  was 
determined  according  to  a  procedure  discussed  in  Beers'  theory 
of  errors  (Beers,  1957) .  The  maximum  measurable  error  in  all 
cases  was  only  a  small  fraction  of  one  percent. 
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Figure  12  A  logarithmic  plot  of  a  typical  data 
record  (solid  line)  .  The  dashed  line  is  a  least 
mean  squares  fit  of  the  data  displayed. 
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Figure  13  A  logarithmic  plot  of  a  typical  data 
record  (solid  line)  .  The  dashed  line  is  a  least 
mean  squares  fit  of  the  data  displayed. 
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Figure  14  A  logarithmic  plot  of  a  typical  data 
record  (solid  line)  .  The  dashed  line  is  a  least 
mean  squares  fit  of  the  data  displayed. 


Pressure  =  500  kPa.  Temperature  Difference  =  398.2  K 

1U'1 

1 

i  102 

1  1 1 1            1 

.3 

!                                     ^^             : 

1 

00 

'                ^^^^^                ' 

Filtered 

\<^          \ 

10-« 

.8 

3             0.1            0.2            0.3            0.4            0.5            0.6            0.7            0 

Time  (sec) 

Figure  15  A  logarithmic  plot  of  a  typical  data 
record  (solid  line)  .  The  dashed  line  is  a  least 
mean  squares  fit  of  the  data  displayed. 
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V.   RESULTS  AND  DISCUSSIONS 

A.   ABOVE  ONSET  DATA 

Data  was  obtained  for  a  helium  filled  prime  mover  above 
onset  for  mean  gas  pressures  of  238,  308,  376,  and  500  kPa, 
and  differential  temperatures  ranging  from  onset  up  to 
approximately  400  K,  in  10  degree  increments.  The  pressures 
were  chosen  because  they  cover  the  range  of  the  " ideal"  gas 
pressure  for  the  plate  spacing  in  the  prime  mover  stack  and 
because  they  match  the  pressures  analyzed  in  Lin's  1989 
thesis.  Although  Lin  also  analyzed  the  pressure  of  170  kPa, 
it  was  not  used  in  this  thesis.  This  is  because  the  onset 
temperature  of  the  prime  mover  at  170  kPa  is  so  high  that  an 
insufficient  range  of  data  could  be  measured  with  a  self 
imposed  hot  end  temperature  maximum  of  7  00  K. 

The  results  for  the  four  pressures  are  presented  in 
Figures  16  through  19  which  show  1/Q  as  a  function  of  the 
temperature  difference  across  the  stack.  As  can  be  seen,  the 
collected  data  closely  follows  the  slope  of  the  theoretical 
curves  and,  for  high  mean  gas  pressure,  the  data  and  theory 
exactly  agree.  As  the  mean  gas  pressure  decreases,  however, 
the  actual  value  at  which  onset  occurs  becomes  increasingly 
higher  than  the  theoretical  value.  Consequently,  the  spread 
between  the  actual  data  and  the  theory  increases  as  well. 


35 


-0  001 

1 

+ 

1                               1 

* 

1 

- 

1 238  kPi| 

-0  002 

_                                \ 

+ 

g 

-0.003 

+ 

+ 

-0.004 

\                                                 + 

1 

1           \               1 

+ 
1 

""^20 

340 

360                        380 

Temperature  Difference  (K) 

400 

420 

Figure  16  Data  obtained  for  238  kPa  plotted  as  a 
function  of  the  temperature  difference  across  the 
stack.   The  solid  line  is  the  theoretical  result. 
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Figure  17  Data  obtained  for  308  kPa  plotted  as  a 
function  of  the  temperature  difference  across  the 
stack.   The  solid  line  is  the  theoretical  result. 
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Figure  18  Data  obtained  for  376  kPa  plotted  as  a 
function  of  the  temperature  difference  across  the 
stack.   The  solid  line  is  the  theoretical  result. 
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Figure  19  Data  obtained  for  500  kPa  plotted  as  a 
function  of  the  temperature  difference  across  the 
stack.   The  solid  line  is  the  theoretical  result. 
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The  results  can  also  be  examined  from  a  different  point  of 
view.  Instead  of  an  absolute  comparison  of  theory  and  data, 
it  is  useful  to  consider  1/Q  as  a  function  of  the  number  of 
degrees  above  onset.  In  this  way,  the  actual  results  and  the 
theory  simply  start  from  the  common  reference  point  of  zero 
quality  factor.  This  will  consequently  eliminate  any 
discrepancies  that  result  from  failure  to  predict  the  precise 
onset  temperature.  Figures  2  0  through  2  3  present  the  data  and 
the  theory  in  this  fashion.  To  determine  the  onset 
temperature  for  the  measured  data,  a  least  mean  squares  fit 
was  performed  on  the  data  and  the  result  was  extrapolated  to 
find  where  1/Q  equaled  zero.  What  these  plots  reveal  is  that 
at  high  mean  gas  pressures,  theory  exactly  agrees  with  the 
measured  data  over  the  entire  range  of  differential 
temperature  measured  (approximately  80  K)  .  Again  however,  as 
the  mean  gas  pressure  decreases,  the  data  visibly  deviates 
from  the  expected  theory  within  this  differential  temperature 
range  observed. 

Summarizing,  the  data  displays  two  effects  relating  to  the 
mean  gas  pressure.  First,  as  the  mean  gas  pressure  decreases, 
the  theoretical  value  of  onset  is  increasing  less  than  the 
measured  value;  and  second,  as  the  mean  gas  pressure 
decreases,  the  slope  of  the  actual  data  is  increasingly  less 
than  the  theoretical  slope.  Three  possible  explanations  exist 
as  to  why  these  deviations  occur.  First,  imperfections  in  the 
construction  of  the  prime  mover  may  cause  the  prime  mover's 
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Figure  20  Data  obtained  for  238  kPa  plotted  as  a 
function  of  the  number  of  degrees  above  onset.  The 
solid  line  is  the  theoretical  result. 


Figure  21  Data  obtained  for  308  kPa  plotted  as  a 
function  of  the  number  of  degrees  above  onset.  The 
solid  line  is  the  theoretical  result. 
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Figure  22  Data  obtained  for  376  kPa  plotted  as  a 
function  of  the  number  of  degrees  above  onset.  The 
solid  line  is  the  theoretical  result. 


0.000  m 

1 

1                              | 

1 

-0.001 

[500  kPa| 

-0  002 

_ 

2 

-0.003 
-0.004 
■0  005 

I 

1                             I 

1 

( 

) 

20 

40                          60 

80 

1( 

X> 

Degrees  Above  Onset  (K) 

Figure  23  Data  obtained  for  500  kPa  plotted  as  a 
function  of  the  number  of  degrees  above  onset.  The 
solid  line  is  the  theoretical  result. 
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output  to  be  less  than  ideal.  Second,  a  difference  in  the 
rate  of  harmonic  growth  could  be  responsible  for  the  data 
which  is  deviate  from  the  theory.  And  third,  the  deviate  data 
could  have  resulted  if  this  data  was  measured  during 
conditions  of  higher  thermoacoustic  heat  transport. 

Since  the  prime  mover  is  not  of  ideal  construction,  its 
performance  would  be  expected  to  be  less  than  the  theoretical 
prediction.  This  is  what  has  been  observed.  In  all  cases, 
the  measurements,  at  best,  matched  the  theoretical  results. 
The  actual  energy  transport  taking  place  in  the  prime  mover 
stack  occurs  in  a  fluid  channel  which  is  approximately  one 
thermal  penetration  depth  thick  (<SJ  and  centered 
approximately  one  thermal  penetration  depth  from  the  plate 
(Swift,  1988,  p. 1152).  This  is  important  because  as  the  mean 
gas  pressure  is  decreased,  6K  increases.  Since  SK  is  also  a 
function  of  temperature,  its  value  changes  along  the  stack. 
For  the  two  lowest  pressures  looked  at  though,  SK  at  the 
average  stack  temperature,  is  greater  than  the  half -gap 
spacing  of  the  plates,  and  for  the  two  highest  pressures,  it 
is  smaller.  Consequently,  any  imperfections  in  the 
arrangement  of  the  plates  would  increasingly  interfere  with 
the  prime  mover's  performance  as  the  mean  gas  pressure  is 
decreased.  At  higher  mean  gas  pressures,  the  misalignment  of 
the  plates  simply  becomes  less  of  a  concern  because  the  heat 
transport  occurs  closer  to  the  plates. 
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The  theory  does  not  take  into  account  energy  lost  to 
harmonics.  So  if  harmonics  are  present,  they  would  result  in 
the  absolute  value  of  1/Q  being  less  that  predicted.  Analysis 
reveals,  however,  that  harmonic  growth  is  not  the  cause  of  the 
deviations.  This  is  evident  because  the  growth  of  higher 
harmonics  was  observed  to  be  relatively  independent  of  the 
mean  pressure.  Also,  harmonic  generation  should  be  stronger 
for  larger  acoustic  amplitudes.  Therefore,  it  would 
reasonably  be  expected  that  the  higher  AT  data  would  be  more 
effected  than  the  lower  AT  data  since  /3  is  much  larger.  It 
can  be  seen  however,  from  Figures  20  and  21,  that  the  slope  of 
the  measured  data,  though  in  deviation  with  the  theory,  is 
very  linear.  No  AT  dependence  is  evident.  Also,  we  were 
careful  to  exclude  data  records  which  had  noticeable  harmonic 
content. 

Lastly,  the  question  of  heat  transport  must  be  addressed. 
The  entire  reason  for  looking  at  signals  of  low  amplitude  is 
to  ensure  that  the  heat  transported  by  the  sound  is  not  yet 
significant.  This  is  desired  because,  after  the  heat 
transport  becomes  significant,  the  heat  exchangers  initially 
have  trouble  keeping  up  with  the  energy  demand  placed  upon 
them.  The  transport  of  heat  through  the  fluid  is  simply 
faster  than  the  transfer  of  energy  through  the  solid  heat 
exchangers.  Consequently,  the  temperature  of  the  hot  and 
ambient  end  begin  to  change  and  the  stack's  condition  becomes 
dynamic.   With  thermal  equilibrium  lost,  the  measured  values 
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no  longer  represent  j3.  Since  heat  transport  is  proportional 
to  the  square  of  the  amplitude  of  the  pressure  wave,  it  is 
important  to  know  at  what  pressure  the  signal  was  at  when  (3 
was  measured.  Figures  24  through  27  show  the  range  of  the 
peak  acoustic  pressures  generated  by  the  prime  mover  in  the 
data  records  used  in  determining  the  values  of  /3.  The 
pressures  are  small  enough  that  the  resultant  heat  transport 
is  insignificant.  Also,  they  are  all  approximately  the  same 
order  of  magnitude  for  the  various  mean  pressures.  As  a 
result,  any  heat  transport  would  be  essentially  the  same  for 
high  and  low  mean  pressures.  Finally,  deviations  due  to  heat 
transport  would  be  AT  dependent  because  the  pressure 
amplitudes  are  higher  for  higher  values  of  AT.  As  noted 
above,  however,  no  AT  dependence  is  evident.  It  can  be 
concluded  therefore,  that  heat  transport  is  not  responsible 
for  the  deviations. 

B.   SOURCES  OF  ERROR 

The  theory  discussed  in  Chapter  II  does  not  easily  handle 
the  insertion  of  the  butterfly  control  valve  in  the  prime 
mover  ambient  end.  To  see  how  much  effect  the  valve  has  on  Q, 
the  Q  of  the  prime  mover  was  measured  with  no  temperature 
difference  applied  to  the  stack.  These  results  were  then 
compared  to  Lin's  measurements  of  Q  for  the  same  prime  mover 
but  without  the  valve.  The  measurement  procedure  used  to 
obtain  Q  was  identical  to  the  procedure  used  by  Lin  in  1989. 
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Figure  24  Range  of  the  peak  acoustic  pressure  of 
each  data  record  used  in  determining  j3.  (P  =  238 
kPa) 
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Figure  25  Range  of  the  peak  acoustic  pressure  for 
each  data  record  used  in  determining  0.  (P  =  3  08 
kPa) 
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Mean  Pressure  =  376  kPa 
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Figure  26  Range  of  the  peak  acoustic  pressures  for 
each  data  record  used  in  determining  /3.  (P  =  376 
kPa) 
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Figura  27  Range  of  the  peak  acoustic  pressures  for 
each  data  record  used  in  determining  /3.  (P  =  500 
kPa) 
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The  results  of  this  comparison  show  that  the  Q's  agree  within 
-1%.  (Lin,  1989,  pp.  21-34) 

The  largest  experimental  error  is  likely  due  to  the 
measurement  in  temperature.  The  error,  based  on 
manufacturer's  specifications,  is  ±2°C  for  type  E 
thermocouples  and  ±5°C  for  type  K  over  the  temperature  range 
of  interest. 

As  discussed  at  the  beginning  of  Chapter  IV,  the  error  in 
mean  pressure  is  ±1  kPa. 
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VI.   CONCLUSIONS 

It  has  been  observed  that  the  theoretical  predictions  of 
the  work  output  of  a  prime  mover  are  in  very  good  agreement 
with  the  data  for  high  mean  gas  pressures.  As  mean  gas 
pressure  decreases,  however,  the  agreement  diminishes. 
Likewise,  the  prediction  of  the  slope  of  the  data  above  onset 
is  very  good  at  high  mean  gas  pressures.  It  is  exact  over 
large  ranges  of  differential  temperature.  Again  however,  as 
the  mean  gas  pressure  decreases,  the  slope  of  the  theoretical 
values  of  1/Q  increasingly  exceed  those  of  the  experimental 
results.  These  results  are  consistent  with  those  found  for 
the  same  prime  mover  below  onset  (Lin,  1989) .  As  explained  in 
the  previous  chapter,  these  observed  deviations  are  probably 
due  to  the  inevitable  construction  irregularities  of  the  prime 
mover  stack.  This  method  for  measuring  1/Q  has  proven  to  be 
very  effective  for  low  amplitude  signals. 

To  verify  the  discussion  of  the  results  in  the  previous 
chapter,  it  is  recommended  that  one  of  two  things  be  done. 
First,  a  working  gas  other  than  helium  could  be  chosen.  As  a 
result,  the  thermal  penetration  depth  6K  would  be  different 
from  those  employed  in  this  thesis.  If  the  same  SK  to  plate 
spacing  ratios  are  examined  as  those  in  this  thesis,  then 
similar  results  are  expected.  Second,  a  new  prime  mover  stack 
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could  be  constructed  with  a  plate  spacing  different  from  that 
of  the  stack  used  here.  When  analyzed  at  pressures  that  give 
the  same  penetration  depth  to  plate  spacing  ratio,  similar 
results  are  expected.  If  either  experiment  produces  results 
that  differ  from  those  seen  here,  then  it  would  indicate  that 
the  source  of  the  deviations  observed  in  this  thesis  is 
something  more  than  just  a  physical  construction  problem. 
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APPENDIX  A 


PRIME  MOVER 

DIMENSIONS 

Prime  Mover  Internal  Length 

Prime  Mover  Internal  Radius 

100.05  cm 

1.91  cm 

COMPONENT  DIMENSIONS  AND  LOCATIONS 


Component 

Ambient  Heat 
Exchanger 

Stack 

Hot  Heat 
Exchanger 

Width  (Ax) 

2 . 18*  cm 

3.50  cm 

0.762  cm 

Plate  half- 
gap  (y0) 

0.052  cm 

0.03  8  5  cm 

0.052  cm 

Plate  half- 
thickness  (1) 

0.02  2  5  cm 

0.014  cm 

0.0225  cm 

Total  plate 
perimeter 

312  cm 

2  06  cm 

15  6  cm 

Location  (x) 

88.11  cm 

9  0.29  cm 

93.79  cm 

*  Includes  0.15  cm  gap  between  twin  heat  exchangers. 


49 


APPENDIX    B 


DATA  RECORD 

SUMMARY 

Pressure (kPa) 

Temperature 
Difference (K) 

Range 

containing 

only  first 

harmonic 

Filtered 
range  used 

in  (3 
measurement 

238 

339.6 

1-4000 

300-3700 

238 

349.2 

1-3000 

300-3000 

238 

359.7 

1-4000 

300-3700 

238 

367.2 

1-4000 

300-3700 

238 

377.8 

1-3000 

300-2500 

238 

388.1 

1-3000 

300-3000 

238 

395.8 

1-2000 

300-2000 

308 

332.5 

1-4000 

300-3700 

308 

342.7 

1-4000 

300-3700 

308 

349.4 

1-4000 

300-3700 

308 

357.6 

1-4000 

300-3700 

308 

366.3 

1-4000 

300-3700 

308 

378.4 

1-4000 

300-3700 

308 

388.1 

1-3500 

300-3000 

308 

397.0 

1-2000 

300-1500 
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DATA  RECORD  SUMMARY  (CONTINUED) 


Pressure (kPa) 

Temperature 
Dif ference(K) 

Range 

containing 

only  first 

harmonic 

Filtered 
range  used 

in  |3 
measurement 

376 

339.3 

1-4000 

300-3700 

376 

349.8 

1-4000 

300-3700 

376 

360.3 

1-4000 

300-3700 

376 

367.8 

1-4000 

300-3700 

376 

378.9 

1-3500 

300-3500 

376 

387.9 

1-3000 

300-2000 

376 

398.1 

1-3000 

300-2000 

500 

361.7 

1-4000 

300-3700 

500 

370.3 

1-4000 

300-3700 

500 

379.6 

1-4000 

300-3700 

500 

388.7 

1-4000 

300-3700 

500 

398.2 

1-4000 

300-3700 
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APPENDIX  C 

EQUATIONS  USED  FOR  DETERMINING 

TEMPERATURE  DEPENDENT  QUANTITIES 

Mean  Temperature  as  a  function  of  x; 

Tm[x]    =  kT0  +  kT1{x  -  xsCack)    +  kT2(x  -  xSCack) 


TO    "    ^TOO  *~T01^  ^ 

kTl     =    CT10    +    CTlAT    +    CT12(M^ 
kT2    =    CT20    +    CT2AT    +    CT22^T> 


2 


cT00  =  292.86 
CT01   =  0.00171124 


cT10  =   -12.616  5 

ctii  =  29  •  0362 
cT12  =  0.00920377 


cT2o  =  335  .  803 
CT21   =   -12.3295 
CT22  =   "0.27  0547 


Mean  Temperature  Gradient: 


VTm[x]   =  kg0   +  kgle-k«(x  -  x'^] 


kg0 

Kg2 

= 

Cg00    +    CgoAT  + 
Cgio    +    CglAT  + 

cg22(*T)2 

Cg00 

= 

9.36469 

Cg01 

= 

18.9501 

Cg02 

= 

0.0000435077 

^gio 

= 

-9  .90268 

cgii 

= 

9 .63328 

Cgl2 

= 

0.0122971 
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Sound    Speed 


Cg20   =   0.  1215  56 
cg21    =   0.0722492 
Cg22   =    -0.0000210737 


a  [21    =  a0  JT7T~0 

aQ  =  1008  m/s 
T0  =  29  3     K 


Angular    Frequency; 
w  [ATI    =  


rm 


stack        ,     , .     r  -^  _i_  A  „  i  ^ stack         ""^scacJc 


+    <«>l[*seacJfc  +  A*sCac*]      + 


\  ■*>       "      X_,-,_j,      "      ^X„rxr-lr' 


-,  r  t<        i  _i  L--scacJc stacK1  ,rm        i 


o),  [x]    =  2/!Tj/2^  +  VTlxj  (x  -  xstacjc)    -  yr^]  -    a^^] 


Phase: 

<|>[x]    =  ^[AH^  -  VTLxj  (x  -  xsCacJt)    -  JT^ 
4),  [AT]    =  20>[Ar]N/T^/(aoVT[x]) 


Miscellaneous: 


x[7l  =  (i>[Ar]  /a  [71 

|i[7l  =  \i0(T/T0)* 

k[71  =  K0(r/T0)P 

p. [21  =  p0(p/p0)  (r0/D 

5K[21  =  ^2*121/ (pml2lcPoG)[A2l) 

5V[21  =  ^iilrj/  Cp.lrJwLArJ ) 

a  [21  =    (6AT]/6AT])2 
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Miscellaneous  (continued) : 

P  =  0.694 

Y  =  5/3 

\i0  =  1.961xl0"6  Pa-sec 

K0  =  0.1524  W/m-K 

p0  =  0.164  kg/m> 

PQ  =  1.01x10s  Pa 

cp  =  5197  J/kg-K 
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